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Abstract. For a given crystal structure, say body-centred-cubic, the many-body Hamiltonian H in which 
nuclear and electron motions are to be treated from the outset on the same footing, has parameters, 
for the elements, which can be classified as (i) atomic mass M, (ii) atomic number Z, characterizing 
the external potential in which electrons move, and (Hi) bcc lattice spacing, or equivalently one can 
utilize atomic volume, Q. Since the thermodynamic quantities can be determined from H, we conclude 
that T c , the superconducting transition temperature, when it is non-zero, may be formally expressed as 

T c = T^ M) {Z,n). One piece of evidence in support is that, in an atomic number vs atomic volume graph, 
the superconducting elements lie in a well defined region. Two other relevant points are that (a) T c is related 
by BCS theory, though not simply, to the Debye temperature, which in turn is calculable from the elastic 
constants Cn, C12, and C44, the atomic weight and the atomic volume, and (b) T c for five bcc transition 
metals is linear in the Cauchy deviation C* = (C12 — C44)/(Ci2 + C44). Finally, via elastic constants, mass 
density and atomic volume, a correlation between C* and the Debye temperature is established for the 
five bcc transition elements. 

PACS. 74.62.-c Transition temperature variations - 74. 70. Ad Metals; alloys and binary compounds 

1 Background and outline regarding the role of strong electron-electron interactions 

in both the high-T c cuprates and heavy Fermion systems. 
Here, our basic philosophy will be to insist that if we 



f^j . We have recently been concerned with both empirical and we re able to solve the many-body Schrodinger equation 

£H ■ theoretical relations between the superconducting tran- for the (considered infinite) superconducting materials, 

O ', sition temperature T c of high-T c cuprates and of heavy then by treating the motion of nuclei and electrons on the 

O ■ Fermion materials [1112115]. The generally complex crystal- sa me footing, plus full inclusion of electron-electron inter- 

j>! \ lographic structure of such compounds makes it difficult actions, such uncertainties involved in separating electron- 

*-j . to identify useful correlations between their superconduct- lattice and Coulomb repulsions between electrons would 

' ing properties (such as T c ) and the elastic properties of be bypassed. 
^ \ the lattice. This is not the case of several superconduct- Having said that, let us take as the simplest start- 
ing elements with a definite and relatively simple crys- i ng point the metallic elements. Then, the input infor- 
tallographic structure, e.g. characterized by only a few mat ion into any computer programme to treat these ele- 
non-zero components of the elastic tensor. Although any ments wou ld be as follows. First, of course, we should need 
such correlation applying to the 'simple' superconducting to specify the structure. To be definite, below we shall 
elements may not be immediately generalized to other un- s i ng i e out the body-centred cubic (bcc) lattice, but every- 
conventional superconductors, they are anyway expected thing that follows would be equally applicable to the more 
to focus on the relevant variables which would be worth- c l 0S ely packed face-centred cubic (fee) structure. Once the 
while studying, both experimentally and theoretically, also structure is specified, one would need to insert the atomic 
in the new classes of superconductors. volume Q (or, of course, essentially equivalently, the lat- 
Following the Bardeen-Cooper-Schrieffer (BCS) the- tice parameter a). Then, the external potential created by 
ory 0] of the metallic elements, firmly rooted in electron- the nuclei must be specified, which requires the atomic 
phonon interaction as the basic mechanism resulting in number Z as further input. Since one has a many-body 
the formation of Cooper pairs, questions have come up Hamiltonian containing both electron and nuclear kinetic 
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energies, one needs also the nuclear mass M. Of course, we 
take as obvious the input additionally of the fundamental 
constants h and e, plus the electronic mass m. 

The conclusion from the many-body Hamiltonian is 
therefore that, for a specified structure which we take to 
be bcc for reasons that will emerge below, the supercon- 
ducting transition temperature T c , given from the many- 
body partition function once the Schrodinger equation has 
been solved depends, apart from the given fundamental 
constants h, e and m, on M , Z and ft, that is 

T C = T^ M \Z,Q). (1) 

Of course, for all other classes of superconductors than the 
metallic elements we have more than one atomic number, 
possibly the next simplest case being the alkali-doped ful- 
lerides (see some brief comments in Section |SJ. 

With this as background, the outline of the paper is as 
follows. Sectional picks out specifically five bcc supercon- 
ducting transition elements, W, Mo, Ta, V and Nb. Two 
more elements, Cr and Fe, have low temperature bcc struc- 
tures but exhibit cooperative magnetism at low temper- 
atures (antiferro- and ferro-magnetism, respectively) and 
are not superconductors at the lowest temperature they 
have yet been subjected to. The five elements listed above 
are considered in the (fl, Z) plane with respect to their 
transition temperatures, the reduced isotope effects being 
taken as evidence that in Eq. (|TJ there is, at most, a weak 
and therefore relatively unimportant dependence of T c on 
nuclear isotopic mass. Since even then T c = T^ M) {Z,Q) 
presents problems in its representation, Section [3] reintro- 
duces a classification of the above five elements in which 
T c is related to the Cauchy discrepancy, i. e. the departure 
of C\2 from C 44 , where these are two of the three elas- 
tic constants (Cn being the other) required to character- 
ize a cubic crystal. Section 01 then returns to an essential 
ingredient of BCS theory, and by using a semiempirical 
approach throws light on the way the Cauchy deviation 
relates to the Debye temperature. Section [5] constitutes 
a summary, plus some proposals for further studies, both 
theoretical and experimental, which should prove fruitful. 
An Appendix considers zero temperature properties, and 
in particular critical field H c (0) and energy gap E g (0), as 
functions of the Cauchy discrepancy. 

2 Dependence of T c on atomic number Z and 
atomic volume Q in bcc transition elements 

In Fig. ri] a plot is made of the positions of the five ele- 
ments in the (I?, Z) plane, the values of T c being attached 
to these coordinates. 

That both Q and Z are important variables in char- 
acterizing T c is immediately apparent. As to the func- 
tional form T C (Q, Z), one can comment that: (i) For con- 
stant atomic volume, T c markedly decreases with increas- 
ing atomic number, (ii) For constant Z, there is plainly 
substantial variation of T c with atomic volume, which is 
proportional to the reciprocal of the concentration. Rel- 
evant to such variation is the pressure dependence of T c 
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Fig. 1. Shows five superconducting bcc transition elements 
in the (O, Z) plane. Actual abscissa is the reciprocal of the 
atomic volume, O' 1 . To each point, the value of the critical 
temperature T c is attached. 

Table 1. Pressure derivatives of T c 0, experimental bulk mod- 
uli UJ, and inferred partial derivatives of T c with respect to fi 
at constant Z, Eq. (|5J, of the bcc superconducting transition 
metals at P = 0. 



Element 


W Mo Ta V Nb 


f£ [K/GPa] 
B [GPa] 
§[10 31 K.m-3] 


-1.4 -2.6 6.3 -2.0 
323 272 200 162 170 
2.4 2.9 -7.4 1.9 



for a given element, provided one remains within the bcc 
phase. 

Despite high pressure can turn many elements into su- 
perconductors via an insulator-metal transition, T c usu- 
ally decreases with increasing pressure for most super- 
conducting elements at ambient pressure (see Table 
Within BCS theory [see also Eq. © below] or its ex- 
tension by McMillan, this is usually justified in terms of 
a pressure-induced lattice stiffening, which reduces the 
electron-phonon constant at a more rapid rate than the 
electron density of states at the Fermi level is increased . 
Pressure derivatives of T c can then be straightforwardly 
related to volume derivatives (at constant Z) from the 
relation 

cHogT c _ dlogT e 

dlogfl dP ' U 

where B is the bulk modulus (see Table QJ. 

However, even given some knowledge of these partial 
derivatives, the fact that T c depends apparently in a sen- 
sitive way on these two variables for the chosen bcc struc- 
ture leaves open the detailed form of the function T c {f2, Z) 
for this structure. Therefore in the following section we ap- 
peal to a known, but so far rather neglected, correlation 
between T c and the Cauchy discrepancy Cn — C44 between 
elastic constants. This is important for our present study, 
since it is clear that T c can, in fact, be characterized by a 
single variable, rather than the pair (f2, Z) used in Fig.f^ 
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Fig. 2. Superconducting critical temperatures for bcc tran- 
sition elements correlated with Cauchy discrepancy, Eq. Q. 
Solid line is best fit to the points [see text for discussion, in 
particular Eq. Q]. Redrawn from Ref. |Hj (see also [Sj for an 
update on elastic constants data). 



3 Characterization of T c by the Cauchy 
discrepancy for the five bcc transition 
elements 

Fig. El redrawn from the work of Ledbetter [Hj carried out 
almost a quarter of a century ago, shows a plot of T c versus 
the quantity C* defined by 



C* = 



C\2 — C44 

C12 + a 



44 



(3) 



Ledbetter [Hj also included some alloys, namely Nbo.gZro.i, 
Nbo.4Tio.6, and T10.7V0.3, but we have omitted these from 
the redrawn Fig. even though the alloys support the 
general trend of the correlation shown. Also, the points for 
the bcc elements Cr and Fe have been omitted, since these 
elements are both characterized by magnetic order and 
no superconductivity in normal conditions. It should be 
mentioned, however, that a high-pressure, non-magnetic, 
but also a non-bec phase of iron has been recently reported 
to display superconductivity with T c < 2 K between 15 
and 30 GPa 0. 

The equation of the straight line drawn in Fig. [3 is 



T c [K] = AC* - B, 



(4) 



with A = 17.7 K and B = 1.65 K. Though presently we 
do not have theory to allow the evaluation of A and B 
from first principles, the correlation in Eq. Q leads us, in 
the following section, to attempt to relate C* to a basic 
ingredient of BCS theory, the Debye temperature On ■ 



4 Cauchy discrepancy related to Ob, which 
gives the 'scale' of T c in the BCS theory 

As Allen and Mitrovic ^D] have stressed, notwithstand- 
ing the numerous impressive and successful predictions 
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Fig. 3. Top: Superconducting transition temperatures ver- 
sus Debye temperatures for numerous elements. Redrawn after 
Ref. |11|. Solid and dashed lines are guides to the eye. See text 
for discussion. Bottom: Same as Fig.0 but now with attached 
values of square root of the 'average' elastic constant, as 
defined by Eq. Q, in GPa^. 



for the metallic elements of BCS theory, their formula (see 
Ref. QUI, Eq. (2.29)) 



T c = 1.136> D exp - 



1 



N(E F )V 



(5) 



where N(Ef) is the density of states at the Fermi level 
and V is the electron-phonon coupling constant, is not a 
successful way of correlating values of T c for the metallic 
elements. Nevertheless, it suggests that one should re-open 
empirically the question of a correlation between T c and 
<9d- Therefore, more generally than for the bcc transition 
metals, we have redrawn data by de Launay and Dolecek 
[TT] in Fig. |2| (top panel), adding values for W and Mo. 

While, for mainly non-transition elements, the contin- 
uous line drawn in Fig. [3] (top panel), already given by 
de Launay and Dolecek a decade before BCS theory, 
shows a relation between T c and 6>d, it is far from simple. 
And the triangle involving Ta, V and Nb modified from 
the 1947 figure of de Launay and Dolecek ^I] shows no 
relation to the continuous curve. 

Nevertheless, 6>d lies deeply enough in first principles 
theory to enquire whether it can be connected, albeit not 
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Fig. 4. Shows 'average' elastic constant C, Eq. 0, as a func- 
tion of Cauchy discrepancy C* , Eq. (0 , for the five bcc tran- 
sition metals. Dashed line is a guide for the eye, which ac- 
tually employs an exponential functional form, of the kind 
C = a + /3exp(-7C*), with a = 80.5 GPa, /3 = 1008.9 GPa, 
7 = 12.5. 



simply, with the Cauchy deviation C* , which is much more 
directly related to T c , as shown in the previous section. 

To attempt this, we note that numerous earlier workers 
have calculated the Debye temperature for cubic crystals 
from knowledge of the elastic constants Cn, Cn and C44, 
the mass density and the atomic volume fl. While Hous- 
ton's method fT^] is favoured, and has been developed by 
Bctts et al. |13II14| . we have found the semi-empirical re- 
lation of Blackman |15| . quoted in Huntington's review 
article |16) . a useful starting point. This reads 



u D — 



3.15 

x(C 



11 



c 



12 



/52 Q 

)5(C 11 +C 12 + 2C 44 )^j 



(6) 



where s is the number of atoms in the unit cell and p is the 
mass density. This approximate result, Eq. QBJ, motivates 
the definition of an 'average' elastic constant 



C 



3.15 



pf2: 



S3 



O 2 



(7) 



and Fig. |21 (bottom panel) parallels Fig.^except that co- 
ordinates in the (f2, Z) plane are now labelled by O. Ev- 
idently from this figure at constant volume (1, related 
to 6*d through Eq. J7J increases with increasing Z, in con- 
trast to the behaviour of T c in Fig. ^ Also at constant Z, 
increases with decreasing atomic volume. Neverthe- 
less, again prompted by the BCS theory, we have sought 
to correlate C= with T Cl but now via the Cauchy discrep- 
ancy. Fig. fj] shows, for the five bcc elements, that there is 
indeed a marked correlation, the functional form obtained 
empirically being recorded in the caption. 



5 Summary and directions for future work 

Though the top panel of Fig. |31 makes it quite clear that 
there is no simple relation between Debye temperature 6>d 
and the superconducting temperature T c , we have been 
led, via the considerations from Fig. ri] and Fig. (lower 
panel), to attempt to correlate C, having dimensions of 
an elastic constant and defined in Eq. (JZJ), which in turn 
involved 6>d, with the Cauchy discrepancy C* in Eq. 10. 
These quantities, for the five bcc elements we have fo- 
cussed on here, are clearly inter-related, as Fig. 01 demon- 
strates, and the functional form has been extracted. Since, 
as Ledbetter jS] already pointed out in 1980, T c relates 
linearly to C* as in Fig. there is a clear correlation 
between T c and 6*d, with mass density and atomic vol- 
ume entering through the definition of the 'average' elastic 
constant in Eq. l(7jl. Furthermore, and again motivated by 
BCS theory, zero temperature quantities, namely critical 
field H c (0) and energy gap E g (0) are shown also to cor- 
relate simply with the Cauchy discrepancy C* for the five 
bcc superconducting transition elements in Figs. |S] and El 
The present work stimulates thoughts concerning gen- 
eralization of the basic approach set out here to other 
groups of superconductors. Our view is that the next sim- 
plest class to study is the alkali-doped Ceo compounds, 
the fullcrides, which have been reviewed by Gunnarsson 
|17j . In Fig. 3 of this review, Gunnarsson has plotted T c 
for RbsCeo and K3C60, as a function of lattice parameter, 
which was varied by applying pressure, P. There is a re- 
markably linear increase in T c with P. For Na2Rb 2 ,Csi_ a; C6o 
there is again a linear variation of T c with lattice param- 
eter, but with a much steeper slope. It will be of interest 
for the future to attempt a generalization of the approach 
given here for the metallic elements to the fullerides, in 
view of these striking correlations between T c and lattice 
parameter. 

G.G.N. A. thanks the Department of Physics, University of 
Antwerp, for much hospitality. 



A Zero temperature properties, related by 
BCS theory to T c , as functions of Cauchy 
discrepancy C* 

The purpose of this Appendix is to display a marked cor- 
relation between experimentally estimated values of the 
critical field H c (0) and the energy E g (0), as extracted 
from tunnelling experiments, and the Cauchy discrepancy 
C*, Eq. (j3J. Thus, the upper panel of Fig. shows H c (0) 
plotted against C* for the five bcc transition elements on 
which attention was focussed in the body of the text. The 
data for H c (0) have been taken from Kittel's book JHj (see 
also jT^ ) ■ Evidently, as for T c in Fig. |2 a marked corre- 
lation again exists, the dashed line having the equation 



H c (0) = hC* -b 2 , 
with 61 = 0.41 Tesla and b 2 = 0.05 Tesla. 



(8) 
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Fig. 5. Experimental results for critical field H c (0) and en- 
ergy gap E g (0), as extracted from tunnelling measurements, 
vs Cauchy discrepancy C , as defined in Eq. 
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Fig. 6. BCS ratio E g (0)/kBT c for bcc transition elements. 
Solid line is ratio of Eqs. gj and while dashed line is BCS 
theoretical value. 
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The energy gap E g (0), as estimated from tunnelling 
experiments, has again been taken from data given by 
Kittel ^Hj (see also 0), and is shown in the lower panel 
of Fig. [S] for the same set of bcc transition elements. The 
dashed line, but now with more substantial scatter of the 
experimental points than for either T c or H c (0), has the 
equation 

E g (0) = e x C* - e 2 , (9) 

with e x = 55.8 • 1CT 4 eV and e 2 = 5.16 • 1CT 4 eV. 

BCS theory predicts a constant value for the ratio 
E g (0)/kBT c . Since Eq. |@J and ©, respectively, give a 
reasonable overall fit for the five bcc transition elements 
focussed on in the present study, we have finally plotted 
the 'average' ratio E g {Q)/k^,T c from Eqs. 10} and © as a 
function of C* in Fig. fusing the given parameters A, B, 
ei and e%. 
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